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Abstract 

Motivated in part by string theory, we consider the idea tliat the standard ACDM cos- 
mological model might be modified by the effect of a long-range scalar dark matter in- 
teraction. The variant of this widely-discussed notion considered here is suggested by the 
Brandenberger-Vafa picture for why we perceive three spatial dimensions. In this picture 
there may be at least two species of dark matter particles, with scalar "charges" such that 
the scalar interaction attracts particles with like sign and repels unlike signs. The net charge 
vanishes. Under this condition the evolution of the mass distribution in linear perturbation 
theory is the same as in the ACDM cosmology, and both models therefore can equally well 
pass the available cosmological tests. The physics can be very different on small scales, 
however: if the scalar interaction has the strength suggested by simple versions of the string 
scenario, nonlinear mass concentrations are unstable against separation into charged halos 
with properties unlike the standard model prediction and possibly of observational interest. 
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1 Introduction 



In this paper and its companion |2] we investigate models for a long-range scalar-mediated 
interaction among dark matter particles. Our starting assumptions are that the scalar in 
question is massless, or very nearly so, and that the scalar has negligible effect on physics 
in the visible sector — apart from what results from the gravitational interaction with the 
dark sector. The fundamental theory that motivates our models offers reasons to doubt 
both assumptions, but the arguments are not convincing enough to discourage investigation 
of the possible observational significance of a string-inspired scalar field modification of the 
cold dark matter (CDM) cosmological model. 

In addition to ideas from superstring theory, we are motivated by the thought that there 
surely is considerable room in the dark sector of cosmology for richer physics than the 
standard CDM model. Many current notions about viable examples of richer physics trace 
back to the long history of ideas about scalar-tensor gravity physics, elements of which are 
reviewed in a paper by G. Farrar and one of us The convergence of these long-standing 
ideas about gravity physics with what superstring theory might be suggesting is appealing 
and just might signify that we are being led to a better approximation to cosmology. 

The inspiration from string theory is the Brandenberger-Vafa P picture for why we 
perceive three spatial dimensions. It suggests that the long-range interaction of dark matter 
is the sum of the gravitational force and the force from the exchange of massless scalars and 
gauge bosons. The dark matter particles may be assigned "charges" — to be defined more 
precisely in section |2] — as in electrostatics, but particles with like signs of the scalar charge 
are attracted and unlike scalar signs are repelled. The model requires zero net scalar charge, 
which implies that in linear perturbation theory the evolution of the mass disttribution is 
unaffected by the scalar interaction. This means that initial conditions can be chosen so 
our model and the ACDM cosmology equally well pass the cosmological tests, including 
the new tests based on the measurements of the large-scale distributions of mass and the 
3 K thermal background radiation The scalar interaction introduces a new mode of 
instability, however, in which particles with different signs of charge tend to concentrate in 
different places. If there were a significant primeval amplitude in this charge separation mode 
it would tend to produce observable mass concentrations that are not closely correlated with 
the large-scale mass distribution, an observationally unacceptable situation. This is avoided 
if the initial conditions are very close to adiabatic, as we will assume. If the scalar force were 
comparable to the prediction of the simplest of string theory setups, then even with purely 
adiabatic initial conditions a strongly nonlinear dark matter mass concentration would be 
unstable against separating into massive charged halos, and halos of different signs would 
tend to avoid each other. We will suggest that such curious behavior may have something to 
do with the confusing situation in the theory and phenomenology of structure formation on 
the scale of galaxies. Other viable models for a modified CDM cosmology would have different 

^We take this opportunity to remedy an omission: the survey should have included the work on this 
subject by Frieman and Gradwohl 0]. 
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implications for small-scale structure, of course. An example based on a supersymmetric 
picture for the dark matter, which allows a non-zero net scalar charge, is discussed in our 
companion paper j2]. 

We organize this paper as follows. In section |21 we introduce the Brandenberger-Vafa 
scenario and describe its late-time dynamics in terms of an effective lagrangian for the dark 
matter.^ In section |21 we characterize the long-range forces in the dark sector by a matrix of 
dimensionless quantities which encode the physics of subsequent interest at scales comparable 
to the sizes of galaxies. Using this matrix, we demonstrate in examples how the scalar force 
leads to a charge separation instability in the mass distribution. In section ^ we discuss 
possible observational consequences of this effect for structure formation on the scale of 
galaxies, and we comment on some issues of small-scale extragalactic structure that might 
point to physics along the lines of this model. These considerations are continued in our 
companion paper [2]. In the last section, we return to the fundamental issues mentioned in 
our introductory paragraph: we consider why the scalar force does not manifest itself in the 
visible sector, in the Eotvos experiment and measurements of orbits in the Solar System; we 
discuss the naturalness of assuming such a small mass for the scalar; and we outline some 
other issues of the fundamental physics of the model. 

2 The Brandenberger-Vafa scenario 

Brandenberger and Vafa pQ put forward a simple and intriguing rationale, based on string 
theory, for space to appear three-dimensional. Suppose, the argument goes, the universe 
started off very small and very hot, with all of the nine spatial dimensions predicted by 
string theory curled up at the string scale (which is close to the Planck scale), and all 
the degrees of freedom of string theory thermally excited, including strings wrapping all 
possible directions and strings with momenta in all possible dimensions. If three spatial 
dimensions start to grow large, then the strings wrapping these dimensions will tend to 
annihilate efficiently with one another, because two extended strings moving in three spatial 
dimensions will generically intersect at some time. But if four or more dimensions start 
to grow large, the strings wrapping them will not generically intersect, so they will not 
annihilate one another efficiently, and their tension will draw the size of the dimensions back 
toward the initial state. At late times, then, it makes sense for three dimensions to have 
grown large, the strings extended in these dimensions to have largely annihilated, and the 
remaining six dimensions to have stayed small. 

The strings wrapping and winding the six extra dimensions may stabilize the size of those 
dimensions: winding strings are lighter when the extra dimensions are smaller, but momen- 
tum modes are lighter when the extra dimensions are larger; thus there is a balance point 
depending on the relative numbers of winding and momentum modes. A recent numerical 

^It should perhaps be emphasized that, for us, late time means the matter- and radiation-dominated 
epochs: this makes our investigation somewhat different in spirit from others in the brane gas literature. 
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study of this was presented in 0, for perturbations with no spatial dependence. The basic 
idea of stabihzing moduh of string theory with winding states actually predates P: the 
earliest reference we know is j7j- 

The starting point for the action for dark matter particles in the late-time picture inspired 
by the Brandenberger-Vafa scenario (as in [B]), for many discussions of scalar-tensor gravity 
models, and for our investigation, is 

S = J d'x^g - \{d<pf - V{<P)^^ - Y: J^^ dsm^{<P) . (1) 

The sum is over all particles in the dark sector. For simplicity we ignore the visible sector of 
baryons, electromagnetic radiation and neutrinos; the focus of our discussion is the growth of 
cosmic structure at relatively low redshifts, when the dark matter dominates. We eventually 
take the potential V{(j)) to be a constant, but it is instructive to keep the general function for 
now. For simplicity we will also ignore the cosmological constant (or dark energy), that is, 
we eventually set V{(j)) = 0. At late time in the Brandenberger-Vafa scenario, the winding 
and momentum modes have shed any excitations they carried, and they are heavy (on order 
the string scale in conventional setups) and their interactions are feeble, so it does not matter 
whether they are fermions or bosons, and the world-line description used in (Q) is appropriate 

Neglecting spacetime curvature fluctuations, the wave equation for the scalar field is 

□ — dV/ d(f) — nqdrriq/ dcf) = , ^^2) 

where Ug is the number density of particles of species q and we are assuming the particle 
peculiar velocities are nonrelativistic. The time-scale for relaxation of the value of to the 
minimum of its potential may be comparable to the Hubble time. If (p is still rolling to its 
minimum it can cause serious problems with the cosmological tests jSj. We assume therefore 
that (j) has relaxed close to its minimum or is placed there as an initial condition. 

On scales small compared to the Hubble length, and ignoring V{(j)), one sees from (j2I) 
that the force between two dark matter particles is 

^ Gmim2 ^ drrii dm2 1 
d(f) dcj) A-Kr"^ 

The second term comes from the scalar exchange, computed in linear perturbation theory. 
The drUq/dcj) thus are evaluated at the mean field value. Because of the analogy with the 
Coulomb force law, we will term Qq = dmq/d(f) the scalar charge of particle q. Scalar forces 
are different from the gauge forces familiar from electrostatics: particles with like scalar 
charge sign attract and pairs with opposite sign are repelled.^ Signs are chosen so that F in 

•^It is worth bearing in mind that scalar charge as we have defined it is not actually a conserved charge. 
For example, a heavy particle and its anti-particle have the same scalar charge because their mass is the 
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(jni) is positive when the total force is attractive. 

As we will describe more detail in section ESI the scalar charge-to-mass ratio Qq/mq is 
expected to be on the order of 1/Mpi, where by definition Mp^ = SnG. So the scalar force 
between particles can be expected to be comparable to the gravitational term, independent 
of separation. We see from however, that when V{(j)) is negligible and has relaxed to 
the minimum of its energy the mean scalar charge density vanishes. This means that there 
is no long-range scalar force when the space distributions of the dark matter particles are all 
the same. Thus we can construct viable models in which the scalar force affects the net mass 
distribution only where nonlinear dynamics has driven separation of the scalar charges. 



3 Density perturbations 

We consider here in linear perturbation theory the behavior of inhomogeneous distributions 
of the dark matter components in an expanding universe. The first subsection generalizes the 
usual theory for pressureless freely moving matter to take account of scalar and gauge charge 
interactions among the dark matter components. In section IH.2I we consider the nature of 
solutions to the perturbation equations. The string theory considerations in section lead 
us to our working assumption, as the simplest possibility, that there are two dark matter 
components with masses close to the four-dimensional Planck mass, scalar charges near ±1, 
and gauge charges at a similar scale which however do not substantially affect the dynamics. 



3.1 Linear Perturbation description of the Scalar Interaction 

One can start from the action in equation (^, with a specified functional dependence of 
the masses of the winding and momentum modes on the scalar field value, introduce a 
hydrodynamic approximation, and work out the full set of controlling perturbation equations. 
Anticipating that some readers may be unfamiliar with density perturbation calculations in 
cosmology, we do this in Appendix A, arriving at a special case of equation (jl(J|) below. 
The derivation presented here follows by inspection of the usual analysis of the gravitational 
growth of structure, which we need only modify by adding to the gravitational acceleration 
the acceleration caused by the scalar field inverse square force law. 
We will write the net force between particles q and p as 

Gmprriq _ i , QpQi 

for different species of dark matter particles. As in equation Q, the scalar charge is Qg = 
dmq/d(j), and the condition that the field is locked by its effective potential Vm = ^Uqimq is 

same for every value of 0, so total scalar charge can change via pair annihilation. And if particles with scalar 
charge form a black hole, the scalar charge is lost. 
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that the total scalar charge vanishes, 



where riq is the mean number density. These equations give the full nonlinear description 
of the dynamics of structure formation on scales small compared to the Hubble length (and 
assuming the dark matter particle velocities are nonrelativistic, which is realistic well away 
from black holes, and ignoring baryonic matter, which a useful first approximation). They 
also give an easy guide to the linear perturbation equations for the evolution of the distri- 
butions of the dark matter components, as follows. 

Recall that when pressure and all other nongravitational forces may be ignored the mass 
density contrast, 5m,(x, t) = 5p(x, satisfies the linear perturbation equation 

Sm + —5m = 4:7lGp6m , (6) 

a 

where p{t) is the total mean mass density, a is the scale factor, and dots denote derivatives 

with respect to proper cosmological time. The second term arises from the decay of peculiar 

velocities caused by the expansion of the universe, and the right hand-side represents the 

driving term due to the acceleration of gravity. If there are several distinct species of cold dark 

matter particles, still with only gravitational interactions, then the evolution of each density 

contrast 6q is driven by the same gravitational acceleration. The perturbation equations thus 

are ^ . 

Sg H Sg = inGpSm = 4:nG ^ nipUpSp . (7) 

a p 

The mean mass density in dark matter of type p is pp = Upnip, and the total density contrast 



IS 

= ^PpSp/p = fp^P^ 



where fp is the mass fraction in dark matter of type p: fp = Upnip/ J^g^q^g- 

The scalar interaction changes ((Tj) by adding to the right-hand side of the equation the 
term representing the acceleration of matter of type q by the scalar interaction with matter 
of type p: 

Sg + —6q = AttG Y npSp {nip + QpQg/ AnGrriq) . (9) 
a p 

The ratio of the two terms in parentheses is defined in equation so we can rewrite this 
expression as 

K + —K = 47rGp pQpfp^p ■ (10) 
a p 

This is the desired result. 

It is an easy exercise to obtain equation ()10|) in a more direct way from the linearized 
expressions for mass conservation and the equation of motion, along with Poisson's equations 



5 



for the gravitational and scalar forces, as in jH]. Yet another approach to this result is 
presented in Appendix A. 

The dark matter components may have gauge charges Q| as well as scalar charges Q^. 
Since peculiar velocities are small, and the gauge charge-to-mass ratio suggested by string 
theory is exceedingly small compared to electromagnetism in the visible sector, the gauge 
interaction simply adds another long-range force that brings the matrix in equation (0]) to 

n _ . QpQq " QfQq /-..n 

AnGmpTTig ' ^^^^ 

where we have summed over s and g to indicate that there may be several different scalar or 
gauge interactions. Under each type of interaction, the net charge must vanish. In the case 
of gauge charges, this follows from Gauss's Law, while for the various scalars, the argument 
again is that the effective potential for the scalars is minimized. Thus 

E fpQl/^P = 0' E fpQl/^P = 0, 

(12) 

P 

where the second line follows from the first applied to 

In the following subsection, we will exhibit solutions to pOj) in the regime where non- 
relativistic matter is the dominant component of the energy density of the universe. Let us 
give here a brief intuitive preview of the results and comment briefly on the evolution of 
perturbations in earlier epochs. 

First, there is always an adiabatic mode whose evolution in the linear regime is exactly 
as it would be in the standard ACDM model. To show this, multiply (fTUj) by fq = mqUq/ p 
and sum over q. Second, there are isocurvature charge separation modes: a gauge charge or 
a scalar charge (or both) develops locally, compensated by an opposite charge elsewhere. In 
the case of gauge charges, such charge separation grow less fast than the adiabatic mode. 
This is because in the absence of gravity, a charged plasma dynamically seeks local charge 
neutrality and leads to Debye screening for the photon. Exactly the opposite is true of scalar 
charges because like charges attract and opposite charges repel: isocurvature modes which 
lead to separation of scalar charges can grow faster than the adiabatic mode. We will see 
that a reasonable suppression of isocurvature modes in the initial spectrum of perturbations 
avoids conflict with observation from the linear regime. In section |3] we will consider the non- 
linear regime, where scalar charge separation modes are important even when isocurvature 
modes are suppressed in the initial spectrum. 

We remind the reader that equation (jiup for the evolution of the density contrasts in the 
dark matter assumes all the components are nonrelativistic and all other kinds of matter 
and energy produce negligible peculiar gravitational accelerations. However, since (jlUp uses 
the local physics of relative accelerations it is valid even in the presence of a cosmological 



6 



constant or of a significant mass density in homogeneously distributed radiation. In the 
presence of such terms the mass density p in equation (|lUj) remains that of the dark matter, 
and when p is significantly smaller than the total source density ptot + Ptot that is driving 
the expansion of the universe, the right hand side of pO|) is not a significant source term 
for the evolution of the 5q. This is the situation at redshifts z ^ z^c^ ^ 3000, where the 
mass density in radiation and neutrinos dominates the dark matter, and the pressure and 
self-gravity of the radiation are important. Since the visible sector is not supposed to be 
significantly affected by the scalar field, the behavior of our model prior to domination of the 
expansion rate by the dark matter is very close to ACDM provided the initial scalar charge 
separation is small. 



3.2 Solutions 

Let us now solve the linear equations (|10p in the approximation of neglecting all sources of 
stress energy except the dark matter. Consider a linear combination of the density contrasts: 

^c = Ecg/A'5g, (13) 
1 

where the Cq are at this point arbitrary constants. The factor of ^JYq is introduced so as 
to extract from (fTn|l an eigenvalue problem for a symmetric real matrix, making the com- 
pleteness properties of the solutions obvious. One sees from ()10p that the linear combination 
satisfies 

Ac + — Ac - AnGp Cp^p/3pqfq5q = . (14) 



If we choose the Cq so that 

"^2, ^P\J~fp(^P<l\J~f'<l ~ ^cCq (15) 

p 

for some number C,c, then ()14|) involves only. This is an eigenvalue problem: if we define 

^pg '\l~fpPpq\J~fq 1 (-'-6) 

then the are the eigenvalues of and the Cq are the corresponding eigenvectors. Since 
is real and symmetric, the eigenvectors form a complete set, and we may completely 
decouple (fTUI) into equations of the form 

Ac + —Ac - 47rGpecAe = . (17) 
a 

For a = (t/to)^^^ (as appropriate for the matter-dominated era), the solutions to ()17|1 are 
Ac ~ t27,±/3^ where 

,,, = zi±^I±^. (18) 
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When the square root is real, we will set 7c = 7c+, for the more rapidly growing mode. 

Under the condition that the total scalar charge vanishes, so equation (|12p applies, one 
solution to the eigenvalue problem always is the adiabatic mode, Cq = ^JYq- This is the mass 
density contrast mode, Ac = 5m- Its eigenvalue is = 1, and the rate of growth of the mass 
density contrast, 5m oc t^^^ in the matter-dominated epoch, follows the standard model. The 
power law indices for the evolution of the isocurvature modes are constrained by the trace 
of the matrix (|TH|l . Since = 1, the sum of the other eigenvalues is 

E ec = Tr(s) - 1 = E - 1 = 2 E Mij, (e Qf - E Qf) K ■ (19) 

ol g \ s 9 / 

If there are two particle species and the gauge force dominates then the second eigenvalue 
^2 is negative, and the real part of (fTH|) is therefore negative so that the mode amplitude 
oscillates or decays. If the scalar force dominates, then there are positive charge separation 
eigenvalues corresponding to growing mode amplitudes. 
These results are illustrated by some examples. 

1) Pure CDM: Here there is only one species, and /? = / = ^ = 7 = l. The growth law 
5m ~ ^^^^ lies at the heart of the successes of the CDM model. 



2) Two species with gauge interactions: Here the total force law between particles 1 and 
2 is 

(20) 



Gmim2 Q1Q2 



where F is taken to be positive if the force is attractive, and negative if it is repulsive. 
Let us assume that the qi are all either -\-q or — g, while the rrii are all equal to m. To 
achieve overall charge neutrality, one must have f = 1/2 for both species. Then 

= 1 (21) 




72± = -(^-l±Vl + 24e2 

If ^2 < 1/24 the power law indices 72± for this mode are imaginary, which corresponds 
to plasma waves in the visible sector. Since ^2 < the real part of 72± is always 
negative, meaning the amplitude of the charge separation is decreasing. The growing 
charge- neutral mode, with 71 = 1, thus dominates at late times. 

3) Two species with scalar interactions: This is the model we will focus on the most. 
The total force law between particles 1 and 2 is given by 0, provided the scalar is 
canonically normalized. Let's assume that the mass of one species is m(0), the mass 
of the other is m(— 0), and </> = in the background solution. One needs / = 1/2 for 
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both species (scalar charge neutrahty) so as not to push away from 0. Thus 
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2M\ 



2 ( d\ogm \ f 1 
dS J l-l 



PI 



71 



d\ogm^ 



72 



(22) 



1 + 24^^ 



4) Four species with gauge and scalar interactions: In an compactification, winding 
strings carry one conserved gauge charge, while momentum strings carry different one — 
the gauge group is U{1) xU (1). The masses of winding and momentum strings depend 
oppositely on the scalar (if one is m{(j)), the other is m(— 0)), but they do not depend on 
the charge. If = in the background, corresponding to the self-dual radius, then the 
winding charge and momentum charge are equal in magnitude (in the sense of giving 
rise to forces of equal magnitude). Ordering the species as (winding, anti-winding, 
momentum, anti-momentum), one arrives at 
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(23) 



where we've put gravitational forces in the first term, scalar forces in the second, 
gauge interactions between winding modes in the third, and gauge interactions between 
momentum modes in the fourth. To achieve charge neutrality for the two f/(l)'s and 
avoid pushing away from 0, one needs / = 1/4 for all species. Thus 
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dlogm 




^3 = ^4 



-Q' 



(24) 



The first mode is adiabatic (all species move together). The second retains charge 
neutrality in both gauge groups, but winding strings cluster in one place while momen- 
tum strings cluster in another. The two remaining modes involve charge separation. 
Clearly, these last two modes are subdominant (and may even be oscillatory). The 
most interesting physics, then, is the same as in case 3 above. 
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To summarize, in linear perturbation theory there always is an adiabatic mode in which 
the mass density fluctuations evolve in the standard way. Gauge interactions tend to suppress 
the charge separation modes, but the scalar interaction can cause the fractional charge 
separation to evolve faster than the mass density contrast. The dynamical reason for this 
is that gauge interactions tend to be screened by the plasma of charged particles, whereas 
scalars, as our results demonstrate, can be anti-screened, leading to new instabilities. We note 
in section m that if the charge separation mode amplitudes exceeded that of the mass it could 
produce observable mass concentrations at positions that are not correlated with the large- 
scale mass distribution. This effect is not observed, but it is easily avoided by postulating 
that the initial conditions are close to adiabatic, because the charge separation modes have 
been growing only since the expansion rate became matter-dominated at Zeq ~ 3000. For 
example, we will show in section 13.31 that in a particularly simple string theory setup, the 
adiabatic mode grows as t"^^^ while the scalar charge separation mode grows as t. One can 
avoid potential conflicts with observations of large scale structure by letting 6 for the charge 
separation mode be less than a third of S for the adiabatic mode. This corresponds to 
suppressing the charge separation 5 by a factor of ~ 100 at z ~ 3000. 

We have continued using string-inspired terminology in the examples listed above, but 
clearly our setup is a simple generalization of the CDM model: the (3 parameters effectively 
lift the assumption of the universality of free fall in the dark sector, an effect studied in many 
papers [3]. 



3.3 Guidance on Charges and Masses from String Theory 

Let us now consider the simplest string theory setup and actually compute the f3 parameters. 
Reducing the low-energy effective action of string theory from ten dimensions to four on a 
compact manifold K includes the following manipulation (see for example [H]): 



2k^ 



10 



2^10 



d X d yJGioe 



K 



-2<S>i 



i?io + 4(9$io) 



1 



VolK 
2^ 



d^x^G^e-^^ 
d'^Xy/g 



Ra + {d^4 



-diG,,„d^G'^ 



R-^id^.y + ^diGu.d'G'' 



(25) 



-2*10 



9^j 



-<l>4 



Gi'i 



-9s 



2kI = WttG = 2-Ka'gl . 



The string metric Gmn is not rescaled in reducing from ten dimensions to four: the Gij 
are its non-compact components, and the Guv are its compact components. The scale of 
compactification is close to the Planck scale, so it makes sense to let all the fields in ()25|) be 
independent of the directions in K. The term in the action summarizes massless 
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scalar moduli. There are other massless moduh at tree level, coming for instance from the 
anti-symmetric tensor field Bmn- 

Let us assume than in the unperturbed solution, i^' is a square whose constituent 
circles all have circumference 1'n\fa' Let the six coordinates run from to 1'n\fa' ^ so 
that Guv = Suv Now consider perturbing a single diagonal entry: G99 = e'^'^^^^^\ where 
as usual Mp^ = 8nG in four dimensions. The four-dimensional action in Einstein frame 
includes the terms 



S* = y d^Xy^ 



The circle in the 9 direction has radius R = y/a'e'^^^^'^K The lowest mass Kaluza-Klein state 
with momentum in the 9 direction has mass mi = 1/R, and the lowest mass string wrapping 
the 9 direction has mass m2 = R/a'; thus 



[A 



pqi 




(27) 



Scalar charge neutrality gives /i = /2 = 1/2, so ^1 = 1 for the adiabatic mode and ^2 = 2 for 
the charge separation mode.^ The numbers for the Ppq, the fq, and the C,q are all independent 
of the background value of 0: this is a special feature of the exponential dependence of the 
masses on in ()27p. 

A potential difficulty is that in four-dimensional Einstein frame, the effective value of a' 
depends on the four-dimensional dilaton $4. The presence of wrapping and winding strings 
will tend to drive $4 toward more negative values (weak coupling). For the purposes of our 
simple-minded discussion, we might simply assume that $4 is stabilized by some unspecified 
mechanism. Or we could suppose that there are D-branes wrapping extra dimensions to 
stabilize $4. 

The values of the f3pq, and hence of ^2, depend on the normalization of the scalar kinetic 
term. Such normalizations are difficult to compute reliably in string theory because they 
depend on the Kahler potential, a quantity which is not protected from renormalization even 
in the presence of A/" = 1 supersymmetry (it is protected by A/" = 2 supersymmetry). One 
may also object that the model we have used for ()25|1 is excessively naive. Nevertheless, string 
theory provides a clear motivation for the scalar forces to be roughly as strong as gravitational 
forces: these two forces have the same origin in the higher- dimensional description. 

More generally, what is needed to have .^2 on the order of unity is for the masses of dark 



"^T^ is actually not a suitable compactification manifold because it leads to J\f = 8 supersymmetry in four 
dimensions, which does not admit chiral fermions and cannot be broken spontaneously. Working through 
this example, however, gives a good sense of why the f3 parameters are of order unity. 

^If there is an extra species of non-relativistic matter with mass fraction /a and no coupling to the scalar 
field (for example, baryonic matter), then ^2 — > $2(1 — fa) as compared to the situation where this extra 
species is absent. 
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sector particles to have the form 



mg = mqofg{(f)/Mpi) 



(28) 



where fq is a function whose values and derivatives are order unity and otherwise uncon- 
strained. The value of the constant rrigQ is not important. Thus we are not tied to dark 
matter particles with masses near the Planck scale: if there is a mechanism which explains 
the smallness of m^o, then the form (PHj) is natural for regarded as a modulus of some 
family of vacua. ^ 

4 Constraints from astronomy 

A viable alternative to the standard ACDM cosmology must pass a considerable and quite 
demanding suite of cosmological tests To satisfy these tests we follow the standard 
model in assuming the presence of a cosmological constant, A, or a term in the stress-energy 
tensor the acts like it. The values of the dark matter particle masses and scalar charges in 
our model are determined by the dynamics of a scalar field, but the field value has been 
locked — apart from the small spatial gradients driven by the inhomogeneous dark matter 
distribution — at redshifts well beyond those of observational interest. The mean expansion 
history thus is observationally indistinguishable from the standard cosmology. To avoid 
unwanted (because unobserved) effects on the interaction of the distributions of mass and 
radiation at decoupling, we must assume the scalar charge separation is small at decoupling. 
We satisfy this by assuming close to adiabatic initial conditions. Under all these assumptions 
our model passes the cosmological tests based on observations that probe the mean expansion 
history of the universe and the behavior of density fluctuations on scales comparable to the 
present Hubble length. 

The observation that galaxies are useful tracers of mass (see ^U] and references therein) 
requires us to stipulate that the initial conditions are close enough to adiabatic on the scale 
of galaxies that the amplitude of the charge separation mode was small compared to the 
perturbation to the total mass distribution as the galaxies were forming. This is because an 
initially isocurvature mode, with charge density contrasts 6- = —5+ and no perturbation to 
the total mass density, becomes a mass fluctuation as it grows nonlinear, 6- is driven to zero, 
and 5+ grows to values much greater than unity. The result is a mass concentration — usually 
termed a massive halo — compensated by a surrounding hole. 

An example of this effect might be useful. The simple illustration in figure ^ assumes 
spherical symmetry and two types of dark matter with the same mean mass density and 

= (3 = 2, = 0. This means that the force of attraction of like charges is 

twice that of standard physics and the scalar repulsion of unlike charges just cancels the 
attraction of gravity. The initial density contrasts in the growing modes of the two dark 



^To an extent, (|27|l is an example of (|28|l : if the string coupling is very weak, then ^ Mp\. 
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Figure 1: An illustration of the development of a massive halo out of an initially isocurvature 
charge perturbation, that is, an initially homogeneous mass distribution with a small scalar 
charge separation. The system is spherically symmetric. The dashed curves are the densities 
relative to the mean for the two dark matter species, and the solid curve is the total mass 
density. All curves are plotted at the moment when the dense component has just stopped 
expanding in the center. 

matter components are 

N±{<r)/nV{<r) = 1 ± SiCos'^{ivr/2R), (29) 

where A^±(< r) and V{< r) are the particle number and volume within radius r, and n is 
the mean number density. This produces density contrasts 6± = ±5j at the center of the 
system, while at the outer radius R the contained mass is the same as in the background 
cosmological model. The expansion at the outer radius thus follows the background model, 
which is Einstein-de Sitter in this example. The initial value of the perturbation 6i at the 
center is chosen so that after the outer radius of the system has expanded by a factor of 30 
the central part of the denser component has just stopped expanding and is about to start 
collapsing. The figure shows the density runs of the two dark matter components at this 
time, when the central densities are p+/p+ = (37r/4)^ and p^/p~ = 0.442. The solid curve 
shows the total mass density, with a mass concentration at the center and a compensating 
mass deficit near the edge of the system. The lesson is that massive halos can grow out of 
a primeval homogeneous mass distribution with scalar charge density fluctuations through 
the nonlinear development of the charge separation. 

If on the scale of galaxies the primeval charge separation fluctuations extrapolated to the 
present epoch in linear perturbation theory were comparable to the primeval mass fluctu- 
ations it would produce massive halos at positions that are not closely correlated with the 
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large-scale mass distribution. Perhaps the tightest constraint on this effect is the observation 
that no known type of object tends to populate the voids defined by the concentrations of 
normal galaxies with luminosities comparable to the Milky Way Also striking is the 
consistency between the measured anisotropy of the 3 K thermal background radiation and 
what is predicted by initially adiabatic mass density fluctuations that are traced now by the 
galaxies ([S], fOj and references therein). Thus we must postulate that the initial conditions 
are close to adiabatic on scales greater than and comparable to galaxies. 

Under pure adiabatic initial conditions, gravity drives the formation of massive halos that 
are well mixed and fair samples of all the pressureless dark matter components. A well-mixed 
and neutral massive halo can be unstable against separation of the scalar charges, however. 
A particularly simple case, for two charge components, is = 2, (3+- = 0. This means the 
gravitational and scalar forces between oppositely charged components just cancel, while the 
force of attraction of like charges is twice that of gravity. A positive and negative massive 
halo thus would interact with each other only by the gravity of the visible sector mass of 
gas and stars that pulls on both. Massive halos of the two types would tend to form at the 
same place, out of the same primeval mass distribution, but, depending on the relative sizes 
of the isocurvature and adiabatic initial conditions, and on the mass in the visible sector, 
initially overlapping halos may drift apart. 

At larger scalar charge-to-mass ratio, charge separation is driven by a dynamical insta- 
bility, as one sees by the following consideration. Imagine a bound spherical halo of two 
species of dark matter particles with masses m+ and m„ = /im+, with /i < 1. There are 
equal numbers of particles with scalar charges Q_ = —Q+, so the halo has no net charge. 
We will suppose that the two charge distributions behave as two continuous fluids, and that 
the two fluids have close to the same initial distributions so the total scalar charge density is 
close to zero everywhere. Let the characteristic radius of the fluid distributions be r. Under 
slow adjustments of R the energy of the system varies with the radius as 

E = a{l + fi)/r^ -b{l+fiY/r, (30) 

where a and b are constants. The first term on the right-hand side is the kinetic energy. As 
indicated, the rms velocity scales inversely as r under a uniform adiabatic change of radius. 
The second term is the gravitational potential energy, which is proportional to the square of 
the mass and scales inversely as the radius. The minimum of the energy is at 

r = 2a/6(l +/i), E = -b'^{l + fxf /Aa. (31) 

Let us compare this to the energy of the system when one fluid has expanded to a much 
larger radius; the relevant case is where the less massive component has moved to much lower 
density. Now the energy of the system is dominated by that of the more massive component, 
so 

El = a/rl - f3++b/n, (32) 
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when the characteristic radius of the more massive component is ri. The mass that now 
dominates the energy is down by the factor 1 + /i from the first case, and the gravitational 
binding energy is increased by the scalar interaction, which is represented by the factor 
defined in equation (jl)). The energy Ei energy is minimum at 

n = 2a//3++6, El = -b'^(3l^/4a. (33) 

The charge separated state is energetically favored if 

/5++>(l + /.p. (34) 

We do not know an easy way to see whether the well-mixed state in this example has a 
local energy minimum that would slow the growth of charge separation. Our highly simplified 
analysis might be expected to underestimate the instability by ignoring the binding energy 
of the more dispersed dark matter; a similar though even more crude argument suggests 
that a near spherical halo with radially separated scalar charges is energetically disfavored 
against the formation of two separated halos of opposite signs, one less tightly bound than 
the other. Our conclusion is that if is significantly greater than unity, then well-relaxed 
massive halos have a single sign of the dark matter. 

It should be noted that the physics does not set a length scale for this charge separation 
instability: that depends on size of the strongly nonlinear mass concentrations. Thus charge 
separation might be complete within the effective radii (which contain half the starlight) of 
galaxies, charge separation might be in progress on the scale of clusters of galaxies, and it 
likely would not yet be happening on larger scales. In particular, the large-scale galaxy flows 
would be driven by gravity alone, as in the standard model. Also, under hierarchical structure 
formation we would not expect the formation of well-mixed halos of dark matter with ^ 1 
on the scale of present-day galaxies. Rather, charge separation would commence in the first 
generation of bound dark matter halos, and subsequent merging and accretion would favor 
the growth of separately charged massive halos, likely with envelopes of more diffuse dark 
matter with the opposite scalar charge. 

We conclude this section with a few comments on what might learned from numerical 
N-body explorations of how structure would form in this model, and on the possible guidance 
one might find in the phenomenology. We consider the case 2 ^ ^ 4, where the scalar 
force between unshielded charged halos is important but not very much larger than gravity. 

The smoothly distributed dark matter in a cluster of galaxies might be dominated by 
one sign of the scalar charge at the center and the other sign on the outskirts. That would 
not discourage gravitationally-driven merging of clusters, but it might encourage persistent 
separation into two mass concentrations. It would drive the separation of galaxies in a cluster 
that have their own differently charged massive halos, and it would tend to push the stars 
in a galaxy away from the center of the galaxy massive halo, a possibly observable effect. 

A charge-separated massive halo of a galaxy would be about twice as tightly bound as 
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in standard physics, for given dark matter mass and radius. The effect would be observable 
only indirectly, in differences in the motions of satellites that have their own charged dark 
matter halos and the motions of gas and stars that are not bound to satellite halos. 

The nearest galaxy to the Milky Way, the Andromeda Nebula, is falling toward us, 
likely for the first time. It is likely therefore that both of these dominant members of the 
Local Group formed in relative isolation, so each has close to zero net scalar charge, with a 
massive halo of one sign and an envelope of the other. This would mean that the galaxies 
are approaching under the usual gravitational attraction, and the same would be true of 
the irregular galaxies on the outskirts of the Local Group. The galaxy NGC 6822 is about 
500 kpc from the Milky Way, about half the distance of the Andromeda Nebula, and possibly 
in a region where its orbit has been influenced by the scalar force of dark matter that has 
been pushed away from the cores of the two galaxies. We will be watching with interest 
closer analyses of its orbit, as may be possible with new generations of precision astrometry. 

Pairs of galaxies with oppositely charged massive halos would not merge. Pairs of galaxies 
with separations comparable to their effective radii are not uncommon: there is a significant 
number in the rich cluster MS1054, at z = 0.8 ^2]; Hu and Cowie jJS] find a curious number 
of close pairs with close redshifts in a field at 2; ~ 0.4; Struble and Rood ^1] suggest some 
Abell clusters at low redshift have significant numbers of close pairs; binary giant galaxies at 
the centers of clusters of galaxies are not uncommon fSJ^]; and, finally, there are enough 
close pairs of galaxies in the field at low redshift to make the galaxy two-point correlation 
function a good approximation to a power law down to separations r ~ 10 kpc, comparable 
to typical effective radii, meaning the visible parts of the galaxies are close to touching. In the 
standard model a close pair of galaxies soon merges because of the gravitational dynamical 
drag of the overlapping dark halos. The multiple nuclei of brightest cluster galaxies may be 
remnants of recent mergers of clusters, but the idea that the large number of close galaxy 
pairs in MS1054 formed this way is curious because the relative velocities of the galaxies 
are much larger than the velocities within galaxies. It is conceivable also that the orbits of 
bound systems of galaxies in the field decay in just such a way as to preserve the power law 
two-point correlation function, but that seems to be a curious conspiracy. The scalar charge 
model offers a new line of speculation: perhaps close galaxy pairs with the same sign of dark 
matter have long since merged, while pairs with opposite signs of the scalar charge remain, 
loosely bound by the gravitational attraction of their baryons, and long-lived because the 
dominant drag is from the baryonic matter. The concept is curious, but so is the situation. 

We continue this speculative discussion in j2], where the scalar interaction in the dark 
sector has some different features. 
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5 Naturalness, coupling to the visible sector, and the 
weak equivalence principle 



In order to significantly affect galaxy formation, the range of the scalar force needs to be at 
least on the order of a megaparsec. The corresponding mass is unnaturally small: m^0^ is 
a relevant deformation of the lagrangian, and barring some symmetry one would generically 
expect m to be on the order of the Planck scale. It won't do to have be a Goldstone 
boson, because if it were, the non-derivative tree-level couplings we need to the dark matter 
particles would be forbidden. Supersymmetry is another possibility for preventing (p from 
acquiring a mass through loop effects. But supersymmetry, if present at all in Nature, is 
broken at the TeV scale in the visible sector, and based on the conventional understanding of 
supersymmetry breaking, it is hard to see how it would be preserved to the scale of inverse 
megaparsecs in the dark sector. To be more precise, assume that the fundamental scale 
of supersymmetry breaking is 10 TeV, with gauge mediation to the visible sector and only 
gravity mediation to the weak sector. Then based on the conventional understanding of 
supersymmetry breaking, the mass of cf) from a soft supersymmetry breaking term is on the 
same order as the mass of the gravitino — about 0.1 eV. This is still too big by many orders 
of magnitude. We do not have a solution to this apparent fine-tuning problem, which applies 
equally to this paper and to (not to mention the vast additional literature on light scalar 
fields in various contexts). 

We are however disinclined to rule out a model on grounds that naturalness is violated 
only after supersymmetry breaking. After all, the smallness of the cosmological constant is 
precisely such a violation. 

In the following subsections we will argue that before supersymmetry breaking, a light 
scalar which does not roll need not be in severe conflict with field-theoretic notions of natural- 
ness. At first blush this may appear to amount only to the familiar fact that supersymmetric 
theories can have fiat directions. But supergravity corrections to Yukawa couplings lead in- 
evitably to couplings of (f) to the visible sector, and our main purpose is to show that these 
couplings need not spoil standard tests of the equivalence principle, provided the superpo- 
tential and Kahler potential have certain properties. In j2] we will discuss the extent to 
which these properties can be motivated from string theory. A brief preview of the string 
theory motivation will be given in section 15.31 

5.1 Coupling to the visible sector 

The crucial point regarding coupling to the visible sector is that couplings should start at 
dimension six, with a form suggested by supergravity, suppressed by a factor of 1/M|[. 

The dimension six operators in question include a factor of 0^, so must sit close to in 
the background solution — otherwise there is a dimension five component to the visible sector 
couplings which would spoil the argument. This may seem to be a fine-tuning of the state of 
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the universe/ particularly with dark matter species present whose mass is not minimized at 
= 0. There is less difficulty when the scalar force is screened by the effects of an additional 
dark matter species which is massless at = 0, as proposed in Further discussion of this 
point is deferred to |2]. 

If comes from the size of a compact dimension, such as the in the direction 
considered in section then it is not clear that couplings to the visible sector starting at 
0^ can be arranged. We will comment further upon this in section ESI and we will explore 
an alternative geometric origin of in where such coupling is better motivated. For now, 
we will simply assume the minimal couplings of the scalar to the visible sector allowed by 
low-energy supergravity. 

In the framework of d = 4 = 1 supergravity, these minimal couplings are non- 
renormalizable couplings in Yukawa interactions and F-term contributions to the scalar 
potential. It is possible for there to be no such couplings in gauge interactions and D-term 
contributions to the scalar potential. The underlying reason is that the superpotential is 
a section of a line bundle over the Kahler manifold parametrized by the chiral superfields, 
so there are factors of e^^^^^pi in front of each Yukawa couplings, and similarly an overall 
factor of e^/^^pi in the scalar potential; on the other hand, complexified gauge couplings are 
genuinely holomorphic functions of the chiral superfields, so gauge interactions descend from 
the superspace action in a way that does not involve the Kahler potential. We assume that 
the dependence of K on the chiral superfield that includes the scalar is = + . . .: 
this gives rise to a canonically normalized kinetic term for 0. It is also necessary to suppress 
dimension four terms in the superpotential which are linear in and couple to the visible 
sector: these would give rise to dimension five couplings in the lagrangian linear in 0. 

In the framework described in the previous paragraph, nucleon masses will depend on 
the modulus through the mass dependence of the constituent quarks. For instance, for the 
up quark and the proton. 



where naively, Cp = 2e„ + e^. All these e coefficients are expected to be on the order of m„ 
and rud, that is, 5 to 10 MeV. The analogous coefficient for the neutron will be different 
from because the quark content is different. This dependence of nucleon masses on a 
scalar is similar to the situation for variations of the fine-structure constant, except that the 
mass shifts depend quadratically on 0. In the next subsection we will replay some standard 
arguments to set limits on (0). These limits will also tell us the accuracy to which we 
need the Kahler potential and superpotential to satisfy the properties described above which 
lead to The couplings resulting in (jHHjl as well as some of the estimates to follow are 

reminiscent of the discussion of a "Least coupling principle" in |17j . 

^See however |17) for an elegant self-tuning mechanism. 



niu = rriu + e. 



nip = nip + ep 



u 




(35) 
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5.2 The weak equivalence principle 

Assuming the dependence (I35|) . if = in the background, an exchange of two scalars in a 
loop diagram produces an attraction between two nuclei (or more generally between visible 
objects) which is completely negligible in comparison to gravity: 



. - F - ^ - 3:r^ (36) 



and the last expression is small for r ^ l/3Mpi. If, however, (0) ^ 0, then there is a tree- level 
coupling of the nucleon to the scalar: effectively. 



mp = mp + 2ep^^^ , (37) 
and this gives rise to a force which can again be compared with gravity: 

Recall that the e coefficients and hence the scalar forces would differ by factors of order 
unity for protons versus neutrons. Thus the scalar force would be isotope dependent. Tests 
of the universality of free fall put strong constraints on such forces, which roughly amount 
to 7/ < 10~^2. This translates to (0) < 5 x 10~^Mpi. 

Spatial variation of (p could also modify gravitational effects at a measurable level. As- 
sume there is a spatial gradient of (p in the Solar System. The Sun and the Earth, having a 
different fraction of protons to neutrons, will feel the gradient differently, so the scalar will 
make different contributions to their total accelerations. The difference {Aa)s between their 
total accelerations should be very much less than the acceleration of the earth due to the 
sun's gravity, {Aa)g: 



{Aa)s ~ — 
rrir. 



^2 

V 



Mil 



(Aa), = = ^ ^ i . (39) 

^ Rsun 2000 yrs ^ ^ 



Taking {Aa)s/ {Aa)g ^ 10 as a rough guide to experimental bounds [TH], we find that 

^Ho/10\ (40) 



where Hq is the present day Hubble constant: Hq^ ~ 4000 Mpc. If we assume that a 
characteristic scale of variation for the scalar is the size of the galaxy, around 10 kpc, we 
arrive at a variation of across the galaxy on the order 



(A0)galaxy ~ V /2galaxy | V02 1 < 2 X 10-^Mpi . (41) 
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In fact, in our model, scalar charge separation should have developed within the halo of 
the Milky Way, and as a consequence, (0) in the galaxy and V0 are non-zero. To estimate 
how big they should be, note that with all coefficients f3pq of order unity, the acceleration 
of a test dark matter particle at the Earth's location due to the halo should be roughly 
equal to the Earth's average acceleration toward the galactic core, which is f^/-Rgaiaxy with 
V ~ 200km/s. But a dark matter particle's acceleration can be estimated as 

a.^l^|V0HM ^(^. (42) 

This leads to {</))/ Mp\ ~ v^/c^ 5 x 10""^, comfortably within the limits arrived at just after 
dSHl) and in dHj). 

If there is some linear dependence of quark masses on (f) — say a term cie^^/Mpi added 
to the right side of — then one can rework the above estimates to bound ci. Its effects 
are similar to the effects of non-zero (0)/Mpi, and the bounds come out to |ci| ^ 2 x 10"^. 
So for example, a term in the superpotential of the form j^^^YuHuQlUr would have to be 
suppressed by an additional factor of ci to avoid causing problems. 



5.3 Other particle physics considerations 

There is a long list of additional issues of potential interest or concern for the model we have 
presented. Here we will consider only two. 

First, there is a disconnect between the original string theory motivation, where the 
scalar is the modulus controlling the size of an extra dimension, and the optimistic as- 
sumptions about the suppression of couplings of to the visible sector. Conventional string 
compactifications give rise to some dependence of the grand-unified gauge couplings on the 
moduli. 

What we have in mind as a more plausible picture of (p in string theory is a modulus which 
pertains to dynamics which is localized on the compactification manifold. For example, if 
a compactification includes a D3-brane on top of D7-branes, then before supersymmetry 
breaking there is a fiat direction corresponding to moving the D3-brane off the D7-branes. 
If the visible sector is assumed to arise from dynamics localized elsewhere on the compactifi- 
cation manifold, then couplings of the dark sector scalars to the visible sector will indeed be 
suppressed — at least at tree level. In |2| we describe more fully a class of supersymmetric 
field theories which lead to interesting dark sector dynamics that includes scalar-mediated 
forces. We will also revisit in [2] the troublesome issue of supersymmetry breaking. 

An additional subtlety is that dependence of heavy quark masses on the scalar, as in 
indirectly affects the low-energy gauge couplings, because it changes the location of 
the threshold between one behavior of the running couplings and another. We have not 
given careful consideration to the bounds one might place on couplings to the visible sector 
through this indirect effect. 
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Second, cold dark matter coming from massive states of string/M theory cannot be pro- 
duced as a thermal relic. If inflation is to be part of the model, then alternative production 
mechanisms must be considered (see for example [221 references therein). Such mech- 
anisms allow dark matter particles to be much more massive than the weak scale. Such 
particles would not be detected in experiments designed to detect LSP-like WIMPs. Indeed, 
it is difficult to see how one could naturally couple (p at tree level to ordinary WIMPs without 
at the same time coupling it unacceptably to the visible sector. 

Direct detection of LSP-like WIMPs would of course be encouraging to the view that 
the dark matter is comprised of such particles. But until there is a way to experimentally 
estimate the density of such particles, the possibility remains open that they are at most a 
subdominant component of dark matter, and that more massive and/or more exotic objects 
are the dominant component. As long as this possibility exists, it seems worthwhile to 
explore possibilities motivated from string theory or in other ways. 
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Appendix 



The purpose of this appendix is to start with the action 

/ dsmM (43) 

and derive the appropriate special case of (fTO|l as a late time limit, assuming V{(f)) = 0. 
Working through this explicitly will reveal an additional condition on the functions 
mq{(f)) that gives stability at early times, and it will hopefully serve to orient readers un- 
familiar with density perturbation calculations. We will following the conventions of [23 
except for metric signature, which we take to be mostly plus. 

To treat the particles in the hydrodynamic approximation, we introduce the average 
velocities, the number densities and currents, and the stress tensors associated with each 
species: 

Ni = n,Ui = n.mMUp', ■ (44) 

The velocities are subject to the normalization condition [/^ = — 1. The stress tensor for the 
scalar is 

= d'0<p - g^^ Q(90)2 + . (45) 

The following equations dictate the dynamics of the fluids coupled to the scalar field and 
gravity: 

V,Nl = V,r;. = ~n,d,m, 

□ = ^ + E = (T^, - -g.,T, j . 

Here T*-' is of course the sum of the T^^ and T*-'.* 

Consider linear perturbations around a. k = FRW background: 

= + 50 nq = ng{l + 6q). (47) 

We will assume that the background solution has = for all time, so = 50. This can only 
happen if dV/dcj) = and Y^q'^n-qdrUq/ dcj) = at = 0. The background is matter-dominated 
FRW provided we assume that also = at = 0. 

The perturbations ()47|) couple at linear order to perturbations of the metric, but there 
is a simple way of decoupling them: we may consider separately the cases where the total 
matter density perturbation, 5m = Z]g fq^qi vanishes, and the case where 5q = 6m for all 
q. In the latter case (the adiabatic mode), it is straightforward to check that the scalar 

^It takes a bit of thought to convince oneself of the second of these equations {WiT^j ~ —nqdjuiq). A 
quick intuitive check is that for a static homogenous fluid in flat space (with gravity turned off), the j = 
component of this equation states correctly that —pq = —nqd^niq. The equation we have written is the 
unique one which gives the correct result in this case, is generally covariant, and results in a total stress 
tensor T'-* which is conserved when the scalar equation of motion is satisfled. 



s 
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perturbation may be set to 0, and then ()46p reduces to the standard equations of the CDM 
model. Derivations of (jH)), which is an exact consequence of (|46|) for the adiabatic mode, can 
be found in standard references, and we will not repeat the derivation. 

Assuming then that 5m = Yl,q fg^q = 0, we first observe that the perturbation in T*-' 
vanishes at linear order. This may be intuitively obvious, but let us review the argument. 
The scalar contribution to T*-' is second order because T^^ is quadratic in (j). So 

= Y: 6T!^ = Y: [5{nqmq)U'^W^ + nqmq{6qUlW^ + U16W^)] . (43) 

In the background solution, = for all spatial components a and all species q. Pertur- 
bations 6Uq are allowed to be non-zero only for spatial z, so as to preserve the norm of f/*. 
Therefore, the second and third terms in square brackets can be non-zero only when one of i 
and j is and the other is spatial. But we may employ gauge freedom to set ^T*" = 0: this 
is synchronous gauge. The first term in square brackets can be non-zero only for i = j = 0, 
and then it is 



g q \ g 



nq-r^]6<f). (49) 



In the last expression, the sum in parentheses vanishes by assumption, and the other sum 
is proportional to J^q^q^q^ which also vanishes. Thus we can consistently set 6T^^ = 0, as 
claimed. It follows that we may also set 6gij = 0: this is an isocurvature mode. 

We now perturb ()46|1 with the metric held fixed. The conservation equation for A*"* results 

in 

Sq + do.U^=0, (50) 

and da acting on the equation for diT^^ results in 

dm„ 

(51) 



■■ 2a- dlogniq 

5q H Oq = — dada(l) 

a dcp 

where to get the second line we have used (|5U|). For the scalar equation, itself is a first 
order quantity since it is zero in the background, so we have 

where we have used Uq = pfg/rriq in the background solution. At early times, when dada4> 
is negligible on the left hand side of (jSH), the perturbations are stable provided 



m 



d^V , ^ fq d^m 

q ruq dcj) 



pE£^^>0. (53) 



23 



The numerical study [H] applies to precisely such modes whose wavelengths are well outside 
the horizon. 

At late times, when dada(f) dominates the left hand side of (j52j) . one can combine (j51|) 
and to obtain 

■A 2d- ^dlognipdlogmg 

^p + -5.-pE^^ ^/A 

; (54) 

where we have used jSpq = 1 + 2M|j^^^^^^^^^2^^ and J2q fq^q = 0- Recalling that 

P=A'kGp=-'^, (55) 

we see that (jK^ agrees with (|Tn|l for the isocurvature modes. Since we justified (fTn|) for the 
adiabatic mode as well, its derivation as the late time limit for arbitrary perturbations is 
complete. 
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